Tangling clustering instability for small particles in temperature stratified turbulence 
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We study particle clustering in a temperature stratified turbulence with small finite correlation 
time. It is shown that the temperature stratified turbulence strongly increases the degree of com- 
pressibility of particle velocity field. This results in the strong decrease of the threshold for the 
excitation of the tangling clustering instability even for small particles. The tangling clustering in- 
stability in the temperature stratified turbulence is essentially different from the inertial clustering 
instability that occurs in non-stratified isotropic and homogeneous turbulence. While the inertial 
clustering instability is caused by the centrifugal effect of the turbulent eddies, the mechanism of the 
tangling clustering instability is related to the temperature fluctuations generated by the tangling 
of the mean temperature gradient by the velocity fluctuations. Temperature fluctuations produce 
pressure fluctuations and cause particle clustering in regions with increased pressure fluctuations. 
It is shown that the growth rate of the tangling clustering instability is in 

times larger than that of the inertial clustering instability, where Re is the Reynolds number. Ma 
is the Mach number, Iq is the integral turbulence scale and Lt is the characteristic scale of the 
mean temperature variations. It is found that depending on the parameters of the turbulence and 
the mean temperature gradient there is a preferential particle size at which the particle clustering 
due to the tangling clustering instability is more effective. The particle number density inside the 
cluster after the saturation of this instability can be in several orders of magnitude larger than the 
mean particle number density. It is also demonstrated that the evaporation of droplets drastically 
change the tangling clustering instability, e.g., it increases the instability threshold in the droplet 
radius. The tangling clustering instability is of a great importance, e.g., in atmospheric turbulence 
with temperature inversions. 

PACS numbers: 47.27.tb, 47.27.T-, 47.55.Hd 



I. INTRODUCTION 

Formation of spatial inhomogeneities in the number 
density distribution of small inertial particles in a tur- 
bulent flow (also called particle clustering or preferen- 
tial concentration) has attracted considerable attention 
in the past decades. [1-5] The enhanced number density 
of particles inside the cluster may affect the particle inter- 
actions, their dynamics and collisions, e.g. by enhancing 
the relative particle velocity and the rate of particle colli- 
sions and subsequent possible coalescence and/or separa- 
tion of the particles. The dynamics of particle collisions 
is relevant to many phenomena in nature such as the 
raindrop formation and atmospheric aerosols dynamics 
[6-9], as well as to numerous industrial processes involv- 
ing, e.g., spraying of liquid droplets in diesel and rocket 
engines [10]. Formation of clusters with enhanced num- 
ber density of particles may increase the rate of particle 
coUisions and coalescence [6, 8] and subsequently can sig- 
nificantly modify the size and velocity distributions of the 
droplets in the spray, which may, e.g., affect the combus- 
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tor performance. 

Clustering of inertial particles in a turbulent flow is 
caused by the centrifugal effect, which implies that the 
inertial particles are locally accumulated in regions be- 
tween the turbulent eddies. These regions have a low 
vorticity, high strain rate, and maximum fluid pressure. 
Therefore, turbulent vortices act as small centrifuges that 
push heavy particles to the boundary regions between the 
eddies by the inertial forces creating concentration inho- 
mogeneities. This effect is known as the inertia-induced 
particle clustering [11]. The inertial particle clustering 
in a turbulent flow has been studied analytically [12-20], 
numerically [21-28], and experimentally [2, 29-33]. 

There are two types of particle clustering. The first 
type is source clustering related to the source term in 
the equation for fluctuations of the particle number den- 
sity. Most of analytical studies of preferential concentra- 
tion are related to the source inertial clustering. Another 
type of particle clustering is associated with a sponta- 
neous breakdown of their homogeneous spatial distribu- 
tion due to the clustering instability. [15, 17] Although 
nature of the clustering instability is still not fully under- 
stood, the impact of the clustering instability on many 
processes, e.g. mixing, particle transport, etc., can be of 
great importance for different practical applications. 

In the temperature stratified turbulence the parti- 
cle clustering is affected by turbulent thermal diffusion. 
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[34, 35] This phenomenon causes accumulation of the in- 
ertial particles in the vicinity of the mean temperature 
minimum and results in the formation of inhomogeneous 
mean particle number density distributions. Turbulent 
thermal diffusion is a purely collective phenomenon oc- 
curring in temperature stratified turbulence and result- 
ing in the appearance of a non-zero mean effective ve- 
locity of particles in the direction opposite to the mean 
temperature gradient. A competition between the turbu- 
lent thermal diffusion and turbulent diffusion determines 
the conditions for the formation of large-scale particle 
concentrations in the vicinity of the mean temperature 
minimum. The phenomenon of turbulent thermal diffu- 
sion has been studied analytically [34-40], investigated 
by means of direct numerical simulations [41], and de- 
tected in the laboratory experiments in stably and unsta- 
bly stratified turbulent flows [42, 43], and also observed 
in atmospheric turbulence. [40] 

Particle clustering in the temperature stratified turbu- 
lence can be much more effective than the inertial particle 
clustering in isothermal turbulence. [44] The reason for 
this is that the mean temperature gradient in turbulent 
flow is a strong source of the temperature fluctuations 
which are correlated with the fluctuations of fluid ve- 
locity and pressure. The pressure fluctuations increase 
fluctuations of the particles number density and enhance 
the rate of formation of the particle clusters. Moreover, 
tangling of the mean gradient of particle number density 
(formed by the tTirbulent thermal diffusion) generates ad- 
ditional fluctuations of particle concentrations and con- 
tributes to the particle clustering. 

The steady-state regime of the tangling clustering in 
temperatiire stratified turbulence without excitation of 
the tangling clustering instability (i.e., the regime of the 
source tangling clustering) has been recently studied ex- 
perimentally and theoretically in Ref. 44. It was demon- 
strated that in the laboratory stratified turbulence the 
source tangling clustering is much more effective than a 
pure inertial clustering that has been observed in isother- 
mal turbulence. In particular, in the experiments in oscil- 
lating grid isothermal turbulence in air without imposed 
mean temperature gradient, the inertial clustering is very 
weak for solid particles with the diameter « 10 /im and 
Reynolds numbers based on turbulent length scale and 
rms velocity, Re = 250. In the experiments [44] the cor- 
relation function for the inertial clustering in isothermal 
turbulence is much smaller than that for the tangling 
clustering in non-isothermal turbulence. The source; tan- 
gling clustering was studied in Ref. 44 for inertial par- 
ticles with small Stokes numbers and with the material 
density that is much larger than the fluid density. 

The goal of the present paper is to conduct a com- 
prehensive theoretical study of the tangling clustering 
instability in the temperature stratifled turbulence. In 
this paper we show that in the presence of the imposed 
mean temperature gradient the tangling mechanism con- 
siderably enhances the rate of the clustering instability, 
and for small particles the tangling clustering instability 



in the presence of the imposed mean temperature gradi- 
ent may result in the formation of small-scale inhomo- 
geneities with the number density of particles exceeding 
the ambient average particle number density by several 
orders of magnitude. 

The paper is organized as follows. The large-scale ef- 
fects in particle transport in temperature stratifled tur- 
bulence are discussed in Section II. The governing equa- 
tions for analysis of instability are given in Section III. 
Solutions for the tangling clustering instability without 
(Section IV) and with the source term (Section V) are 
analyzed, and the instability growth rate and saturated 
value of the particle number density inside a cluster are 
determined. In Section VI we take into account an effect 
of droplet evaporation on tangling clustering instability. 
Finally, in Section VII we discuss the implications of the 
instability development and conclusions. 

II. PARTICLES IN TEMPERATURE 
STRATIFIED TURBULENCE 

A. Governing equations 

Advection-diffusion equation for the number density 
np{t, x) of inertial particles in a turbulent flow reads [45, 
46]: 

-^ + '^-{npv) = Dm/^np, (1) 

where Z)„i is the coefficient of molecular (Brownian) dif- 
fusion, v{t, x) is the particle velocity field. We use a mean 
field approach in which the particle number density and 
velocity, the fluid temperature, density and pressure are 
decomposed into the mean and fluctuating parts, where 
the fluctuating parts have zero mean values. Averaging 
Eq. (1) over an ensemble of turbulent velocity fields we 
obtain an equation for the mean number density of par- 
ticles: 

dN 

— + V-{NVp + {nv))= Dm AN, (2) 

where N — (up) is the mean particle number density, the 
angular brackets imply ensemble averaging, {n v) is the 
turbulent flux of particles and Vp is the particle mean 
velocity. To obtain a closed mean-field equation one 
needs to determine the turbiilent flux of particles. The 
equation for fluctuations of the particle number density, 
n{t, x) = np{t, x) — N{t, x), then reads: 

On 

— + V-{nv- (nv)) - DmAn = -{vV)N - NV-v. 

(3) 

B. Turbulent thermal diffusion 

Turbulent thermal diffusion results in formation of a 
nonzero gradient of the mean particle number density 
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ViV in temperature stratified turbulence. [34, 35] The 
physical mechanism of turbulent thermal diffusion is as 
follows. For the particles with the material density Pm ^ 
p, their velocity is determined by 



dv 

H 



u — v 



+ 9, 



(4) 



where u is the fluid velocity field, g is the acceleration 
of gravity, Tg = mp/Gnpuap is the Stokes time for the 
small spherical particles of the radius ap and mass rrip, 
v is the kinematic viscosity, p is the mean fluid density. 
For small Stokes numbers, St = Tg/r^ <C f, solution of 
Eq. (4) has the following form (see, e.g., Ref. 11): 



V = U — Ts 



du 

'dt 



+ {u-V)u - g 



+ 0{t!), (5) 



or introducing the dimensionless units, it can be written 
in dimensionless form: 



V = u — St 



du 

'at 



+ 0(St2), (6) 



where distance is measured in the Kolmogorov viscous 
scale £ri = /Re^^^ and the time is measured in the Kol- 
mogorov time scale Tr, = ro/Re^^^. Here Re = iouo/f 
is the fluid Reynolds numbers, uq is the characteristic 
turbulent velocity at the integral scale of turbulent 
motions and tq = (.q/uq. The terms in squared brackets 
in Eq. (5) describe the difference between the local fluid 
velocity and particle velocity arising due to the small but 
finite inertia of the particle. For the turbulent flow with 
low Mach numbers: V-m « —p^^ {u-'V)p ^ 0. The 
equation for V-t) can be easily obtained from the equa- 
tion (5) and the Navier-Stokes equation: 



--{u-V)p 

P 



P 



(7) 



where p are the fluid pressure fluctuations. Due to iner- 
tia, particles inside the turbulent eddies drift out to the 
boundary regions between the eddies. These arc regions 
with small velocity fluctuations and maximum pressure 
fluctuations. Consequently, particles are accumulated in 
the regions with the maximum pressure fluctuations of 
the turbulent fluid, i.e., V-d cx (rs/p)V^p ^ even 
for the incompressible fluid. For large Peclct numbers, 
when the molecular diffusion of particles in Eq. (1) can 
be neglected, we can estimate V- v oc —dup/dt. There- 
fore, incrtial particles are accumulated (i.e., dUp/dt oc 
— (ts/p) > 0) in regions with maximum pressure of 
turbulent fluid, where V^p < 0. Similarly, there is an 
outflow of inertial particles from the regions with the 
minimum pressure of fluid. In case of homogeneous and 
isotropic turbulence a drift from regions with increased 
(decreased) concentration of particles by a turbulent flow 
of fluid is equiprobable in all directions. 



On the contrary, in a temperature stratified turbu- 
lence, the turbulent heat fliix {u 6) docs not vanish. This 
implies that the fluctuations of fluid temperature, 6, and 
velocity are correlated, and, therefore, fluctuations of 
pressure are correlated with the fluctuations of velocity 
due to a non-zero turbulent heat flux, {u9) ^ 0. The 
increased pressure of the surrounding fluid is accompa- 
nied by the particles accumulation, and the direction of 
the mean flux of particles coincides with the direction of 
the heat flux of temperature towards the minimum of the 
mean temperature. [34, 44] 

Equation for the mean number density N of particles 
reads: 

^ + V- [TV (y + Wg) + = Dm^N, (8) 

where Vp = V + Wg is the mean particle velocity, V is 
the mean fluid velocity, Wg = Tgg is the terminal fall 
velocity of particles. Dm = k'BT /Qirpvap is the coeffi- 
cient of molecular (Brownian) diffusion, fcs is the Boltz- 
man constant, T and p are the fluid mean temperature 
and density, respectively. The turbulent flux of particles, 
p{n) _ (jiy^^ includes contributions of turbulent thermal 
diffusion and turbulent diffusion, i.e., 



p(n) ^ yeS JSJ _ n V7V. 



(9) 



Here « £o Uq is the coefScient of turbulent diffusion, 
is the effective pumping velocity caused by the tur- 
bulent thermal diffusion: 



reff 



-t{v {V-v)), 



(10) 



where r is the turbulent correlation time. Equation (10) 

for the effective velocity has been derived using different 
methods in Refs. 34, 36-40. The expression (10) can be 
obtained in a simple way using the dimensional consid- 
eration. Estimating the left hand side of Eq. (3) as 

^ + ^.^nv-{nv))-D„,V\^^, (11) 

we obtain an expression for the turbulent component n 
of the particle nmnber density: 

n « -T V- (Nv) ^ ~t[N{V-v) + {vV)N]. (12) 

Therefore, the turbulent flux of particles i^/"'' = {vi n) is 
given by the following expression: 

p{n) ^ _jvr (t;, (V- v)) - T {viVj)VjN , (13) 

where the first term in the right hand side of Eq. (13) 
determines the turbulent flux of particles due to 
the turbulent thermal diffusion: — N t {vi {'V • v)) = 
Vf^ N, and the second term in the right hand side of 
Eq. (13) describes the contribution of turbulent diffusion: 
T{viVj)VjN = ViN. 

A detailed analysis [34, 44] , using equation of state for 
an ideal gas with adiabatic index Cp/cv (the ratio of spe- 
cific heats) and applying the identity Tg = p Wg Lp/ P, 
yields the effective velocity in the following form: 

(14) 
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where a is 



a = 1 + 



St In (Re) 



Re Ma^ 



(15) 



Here Lp^ = \VP\/P, Ma = uq/cs is the Mach number 
and Cs is the sound speed. For gases and non-inertial 
particles a = 1. A steady-state solution of Eq. (8) for 
V = is given by the following formula: 



N{z) 

Nn 



T{z) 



exp 



Dm + Dj 



dz' 



(16) 



where A^o = N{z — zo) and To = T{z ~ zq) are the mean 
number density of particles and the mean fluid tempera- 
ture, respectively, calculated at the boundary z = Zq. If 
there is a gradient of temperature in a vertical z direction, 
Eq. (16) implies that small particles are accumulated in 
the vicinity of the mean temperature minimum. This 
causes formation of large-scale inhomogeneous distribu- 
tions of the mean particle number density. 



III. GOVERNING EQUATIONS FOR ANALYSIS 
OF INSTABILITY 



Let us study fluctuations of the particle number den- 
sity. The equation for the two-point second-order corre- 
lation function of the particle number density, $(i, R) = 
{n{t, x)n{t, X + R)) was derived using Eq. (3) and path- 
integral approach in a turbulent compressible flow with 
a small flnite correlation time (see Ref. 15): 

^ = [B{R) + 2U^^\R) ■ V + Dij{R)ViVj] $(i, R) 



+I{R), 

where U^^^R) = (1/2) [U{R) - U{-R) 



(17) 



(18) 



D^iR) « 2 / {v, [0, i{t, x\Q)\ V, [t, i{t, X + R\t)\ ) dr, 
Jo 



(19) 



B{R)^2 / {b[0,^{t,x\Q)]b[T,^{t,x + R\T)])dT, 
Jo 

(20) 

/•OO 

Ui{R) « -2 / {vi[0,^{t,x\0)] b[T,${t,x + R\T)])dT. 
Jo 

(21) 

T 

Here b = divv, D^j{R) is the scale-dependent turbu- 
lent diffusion tensor, Sij is the Kroncckcr tensor, I{R) 
is the source of particle number density fluctuations and 
(...) denotes averaging over the statistics of turbulent ve- 
locity field and the Wiener process w{t). The Wiener 
trajectory ^{t,x\s) in the expressions for the turbulent 



diffusion tensor D^j{R) and other transport coefficients 
is defined as follows: 

^{t,x\s)=x— I v[T,^{t,x\T)\ dr — \/2Dmw{t — s), 

J s 

(22) 

where w(t) is the Wiener random process which describes 
the Brownian motion (molecular diffusion) and has the 
following properties: {w{t))w =0, {wi{t + T)wj{t))jj, = 
rSij, and {. . .)w denotes the mathematical expectation 
over the statistics of the Wiener process. The velocity 
Vi[T,^{t,x\T)] describes the Eulerian velocity calculated 
at the Wiener trajectory. 

Equation (17) with I{R) = and for a delta-correlated 
in time random incompressible (6 = 0) velocity field 
was derived by Kraichnan. [47] For a turbulent com- 
pressible flow with a finite correlation time Eq. (17) was 
derived using a stochastic calculus. [15] In particular, 
Wiener path integral representation of the solution of the 
Cauchy problem for Eq. (1), the Feynman-Kac formula 
and Cameron-Martin-Girsanov theorem were used for the 
derivation of Eq. (17). The comprehensive description of 
this method can be found also in Refs. 36, 37, and 48. 

The source function I{R) in Eq. (17) is related to the 
two source terms —{v-'V)N ~ N^-v in the right hand 
side of Eq. (3), and the explicit expression for I{R) is as 
follows (see Ref. 44): 



I{R) = B{R)N'^ + U'^^\R) ■ VAT^ 
+ lD.JR)iS7,N) (V.-AT), 



(23) 



'V)">N{t,x)N{t,y) 



(1/2) [U{R) + Ui-R)] and 
= (3/4) (V,iV) (VjN), and 



where U'^^^ (R) 

r7(a:)Y7(y) 

Vi ^ ^ 

The meaning of the turbulent transport coefficients 
B{R) and U{R) is as follows. The function B{R) is 
determined by the compressibility of the particle veloc- 
ity field. The vector U{R) determines a scale-dependent 
drift velocity which describes transport of fluctuations 
of particle number density from smaller scales to larger 
scales, i.e., in the regions with larger turbulent diffusion. 
The scale-dependent tensor of turbulent diffusion D^j {R) 
is equal to the tensor of the molecular (Brownian) dif- 
fusion in very small scales, while in the vicinity of the 
maximum scale of turbulent motions it coincides with 
the tensor of turbulent diffusion. It should be noticed 
also, that if VA^ ^ a nonzero source term I{R) causes 
the production of the particle number density fluctua- 
tions due to the tangling of the mean particle number 
density by the turbulent velocity field. 



A. Degree of compressibility 

If the turbulent velocity field is not delta-correlated 
in time (e.g., the correlation time is small yet finite). 
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the tensor of turbulent diffusion, D^j{R), is compress- 
ible, i.e., {d/dRi)D^j{R) ^ 0. The parameter that 
characterizes the degree of compressibility of the tensor 
of turbulent diffusion, is defined as follows: 



(7rr = 



^. ((v-l)^) 



(24) 



where Cijk is the fully antisymmetric Levi-Civita unit ten- 
sor, ^ = ^ — x with \t — s\ tq. If the turbulent velocity 
field is a delta-correlated in time random process, then 
V-^ = -(y-v) {t-s), and Vx^ = -(Vx-y) {t-s), and, 
hence, the expression for reads: 



^ ~ ((Vxt;)2) -'^^ 



(25) 



where (Tv is the degree of compressibility of the particle 
velocity field. The parameter (Tv depends on the Stokes 
number [17]: 



(8/3)St 
1 + /3St ' 



(26) 



where the coefficient (3^1. 

For small finite correlation time (i.e., for small Strouhal 
numbers, Sr = Tg-^ (tt^) /£ <^ 1, the parameter can be 
estimated using Eq. (C12) in Ref. 15 (see also Ref. 44) 
as follows: 



The mechanism of coupling related to the tangling of 
the gradient of the mean temperature gradient is quite 
robust. The properties of the tangling are not very sensi- 
tive to the exponent of the energy spectrum of the back- 
ground turbulence. Anisotropy effects do not introduce 
new physics in the clustering process because the main 
contribution to the tangling clustering instability is at the 
Kolmogorov (viscous) scale of turbulent motions, where 
turbulence can be considered as nearly isotropic, while 
anisotropy effects can be essential in the vicinity of the 
maximum scales of the turbulent motions. Using these 
arguments, we consider the tensor D^j{R) for isotropic 
and homogeneous turbulent flow in the following form: 



913 



DiAR) = D, 



+- 



RF' 



Sij - 



F{R) + F,{R)]d,j 
RiRj 



RF, 



f RiRj 



(28) 



where F{Q) = l- Fc(0) and F' = dF/dR. The function 
Fc{R) describes the compressible (potential) component, 
whereas F(R) corresponds to vortical (incompressible) 
part of the turbulent diffusion tensor. 



B. Derivation of expression for the function B{R) 

Taking into account the equation of state of an ideal 

gas wo obtain p/P = g/p + 9/T + 0{gd/pT). For small 
Stokes numbers, V-v {ts/ p)^^P- This allows us to 
estimate B{R) as 

n2 



B{R) 



P 

2Tf P2 
p2 T2 



[VM^)] [V^p(y)]) 
{T[V'e{x)] [V'eiy)]), (29) 



where V2p(a;) = [V^^^] p{x), and we remind that 
the fluid density, temperature, and pressure have mean 
components, p,T,P, and fluctuations, Q,d,p, respec- 
tively. Hereafter we omit the argument t in the cor- 
relation function. In k space the correlation function 
(r[V26l(a;)] [V^0{y)]) reads: 

{T[V'e{x)] [V2%)])= J f{k)kU0{k)e{-k)) 

X exp (ik-R) dk. (30) 

Taking into account that (^(fe) 9{—k)) = 
{6^)E0{k)/A7rk'^, and integrating in k space we ar- 
rive to the following expressions for the functions 

B{Ry. 



BiR) 



2 St^ 4 

3 T^i Uq 



Re 



1/2 



(31) 



where Ee{k) = {2 / "i) k^^ {k / k^) is the spectrum 
function of the temperature fluctuations for ko < k < 
i-^, with fco = ^0 ^ and f(fc) = 2to (fc/fco)"^/^. To deter- 
mine {6"^) we used the budget equation for the tempera- 
ture fluctuations Eg = (0^)/2: 

DE, 



Dt 



-hdiv*e = -{F-V)T-S0, 



(32) 



which for homogeneous turbulence in a steady state 
yields: 

(e^) = -2 TO (F- v)r = I {iovrf , (33) 

where F^ = {uiO) = -D^^^ViT is the turbulent heat 
flux, Dj^^ = is the turbulent diffusion of the tem- 

perature fluctuations and the dissipation rate of Eg is 
eg = (6|2)/2to (see, e.g., Ref. 49). 

On the other hand, the function B{R) at very small 
scales has a universal form: 



B{R) 



2O0-V 



20 cTv 160 St; 



cff 



3t„ 



(34) 



Trj (1 + Tv) uifj 

because at this scales the velocity fleld is smooth and 
nearly isotropic. Here we introduced the effective Stokes 
number using Eqs. (31) and (33): 

Re^/^ 

Steff = St r, r(Ma, Re, Io/Lt) 



9Ma^ 



(35) 
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where Lt is the characteristic scale of the mean temper- 
ature variations. To derive Eq. (34) we took into account 
that for a Gaussian velocity field [15] : 



((V.t,f) = ^St,V 

Equation (36) yields: 

_ _ {{V-vf) 



{{yxvf)=^^. (36) 



((Vxi;)2) 



3 St^ff « 1. 



(37) 



Taking typical parameters for atmospheric turbulence 



Re = 10^ 



Mo 



1 m/s and = 100 m with the tem- 



perature inversion, such that the mean temperature gra- 
dient, I VT|, is of the order of 1 K per 100 m, wc obtain 
r « 2.5 X 10^, or for Re = 10", uq = 0.3 m/s and 4 = 30 
m the parameter F is F w 5 x 10^. 



IV. SOLUTION FOR TANGLING CLUSTERING 
INSTABILITY WITH / = 



It is convenient to rewrite Eq. (17) for the two- 
point second-order correlation function ^{R) in a non- 
dimensional form with coordinate R measured in units 
of the Kolmogorov scale, £ri, time in units of the Kol- 
mogorov time, Tj, and the function $ in units of N^: 



'dt 
where 



1 



$" + 2 ( - 



M(i?) 
+B{R) $ + I{R) 



x{R)] 



2UR^' 

(38) 



M{R) 

xiR) 



S^+3l^ 



F-{RF,y], 



M{R) 



{F - 2F,y, 



I{R) = B{R) + 



4(a- 1)2 



T 

_J_ RxjR) 

2M{R) M{R) 



(39) 



(40) 



Sc = v/Dm is the Schmidt number and U = U{R) R. 
Typically, in many applications, e.g., in the atmospheric 
turbulence, Sc 1 for small inertial particles. The two- 
point correlation function $(/?) satisfies the following 
boimdary conditions: ^'{R = 0) = and $(i? — > oo) = 
0. This function has a global maximum at i? = and 
therefore it satisfies the conditions: 

$"(i? = 0)<0, $(i? = 0) > |$(i?>0)| . 

For a steady-state regime and in the absence of the 
tangling clustering instability the solution of Eq. (38) 
was obtained in Ref. 44 for I{R) ^ 0. In this sec- 
tion we will consider solution of Eq. (38) for the case 
of the tangling clustering instability and I{R) = 0. Since 



the Schmidt number, Sc = f/Dm ^ 1, the molecu- 
lar diffusion scale is much less than the viscous Kol- 
mogorov scale. This implies that the viscous range of 
scales essentially contributes to the tangling clustering 
instability. A general form of the turbulent diffusion 
tensor in the viscous range of scales is obtained tak- 
ing into account that F{R) = (1 — R^)/{1 + a^) and 
Fe{R) = cr^(l - i?^)/(l + cTj.), which yields: 



DiAR) = CiR^5,j + C2R^R, 



2(2 + a^) 

3(1 + C7^) ' 



Co 



2(2(7^ - 1) 

3{l + (7^) 



(41) 
(42) 



and the other functions in this range of scales are U = 

20crv/3(l + (7v) and B = 20a^/{l + a^). 

In the molecular diffusion range of scales, Op/iri < R < 
1/VSc, all terms oc R^ are small and can be neglected. 
Note that we consider the case when the particle radius, 
Up, is the minimum scale of the problem, so that Up < £d, 
where £d = ^o/Pe^^^ = £n/VSc is the molecular diffusion 
scale, and Pe — UQio/Dm is the Peclet number. The 
solution of Eq. (38) in this range reads: 

m)={l-'-^^^^R')eMlt), (43) 

where B > jTr, and 7 is the growth rate of the tangling 
clustering instability. 

In the turbulent diffusion region of scales, 1/VSc ^ 
i? <C 1, the molecular diffusion term oc 1/Sc is negligible, 
and the solution of Eq. (38) in this region we seek for in 
the following form: 



$(i?) cx R-^ exp(7t) 



(44) 



Using the Corrsin integral, R^^{R) dR = 0, we ob- 
tain that /3 = A ± ?'k is a complex number, > 0, where 

^ 4(ij - 7T^)(Ci + C2) - (Ci -C2+ 2Uf ^^^^ 



X = 



4(^1 + (72)2 

C1-C2 + 2U 



2{Ci+C2) 2(l + 3c7^) 



1 + (Tv V 1 + Scr, 



(46) 



Hence the real part of solution (44) is reduced to 

$(i?) = CR-^ cos {KlnR + if) exp(7i), (47) 

where C is the constant. 

Since the correlation function $(i?) has a global max- 
imum in i? = ap/£ri <C 1, the parameter ax < 3. The 
function $(i?) sharply decreases with the increase of R, 
for i? » 1. The growth rate of the second moment of par- 
ticles number density and the constant C can be obtained 
by matching the correlation function ^{R) and its first 
derivative $'(i?) at the boundaries of the above ranges, 
i.e., in the points R = l/\/Sc and R = 1. The matching 
yields k/2(Ci -|- C2) « 7rm/lnSc (where m = 2k + 1), 



FIG. 1. The growth rate 71 of the tangling clustering in- 
stability (in units of l/t,7) of the first mode (m = 1) versus 
the particle radius ap for different values of parameter F, and 
(Ty = 1, Sc = 10® Op. The particle radius Up is given in /jm. 
The dashed line corresponds to the inertial clustering insta- 
bility. 



FIG. 2. The growth rate 71 of the tangling clustering insta- 
bility of the first mode (rn — 1) versus the particle radius Up 
for larger values of parameter F, and a^, = 1, Sc = 10® Op. 
The particle radius Up is given in fim. 

gradient is shown in Fig. 1 by the dashed line). 



and the growth rate for the mode 
clustering instability is given by 



m" of the tangling 



1 



3{1 + 3(T^) 
27r2m2(l 



(H 

3(7^)2 



2(1 + (7^) 



{l + a^)\n Sc 



(48) 



where m = 1,2,3,.... The first mode {m = 1) has the 
minimum threshold for the excitation of the tangling 
clustering instability. 

The tangling clustering instability depends on the ra- 
tio fTy/(Tv. For the 5-correlated in time random Gaus- 
sian compressible velocity field = CTv (for details, see 
Rcfs. 15). In this case the second moment $(t, i?) can 
only decay, in spite of the compressibility of the veloc- 
ity field. On the contrary, for the finite correlation time 
of the turbulent velocity field 7^ Cv, and the correla- 
tion function $(i?) grows exponentially in time, i.e., the 
tangling clustering instability is excited. 

Figures 1-2 show the growth rate (48) of the tangling 
clustering instability of the first mode (m = 1) versus 
the particle radius ap for different values of parameter 
r. One can see from these figures that the characteristic 
time of the tangling clustering instability is of the or- 
der of the Kolmogorov time scale (the growth rate 71 in 
Figs. 1-2 is measured in units of the inverse Kolmogorov 
time). Remarkably that for every parameter r(Re, Ma) 
there is a rather sharp maximum of the function 71 (a^). 
This implies that depending on the parameters of the tur- 
bulence there is a preferential particle size for which the 
particle clustering due to the tangling chistciring insta- 
bility is much faster than for other values of the particle 
size. Moreover, the growth rate of the tangling clustering 
instability is much larger than that of the inertial cluster- 
ing instability (the growth rate for the inertial clustering 
instability in turbulence with a zero mean temperature 



V. SOLUTION FOR TANGLING CLUSTERING 
INSTABILITY WITH 7 

In this Section we obtain solution of Eq. (38) which 
includes both, the tangling clustering instability and the 
source term for the tangling clustering. This implies that 
we consider solution of Eq. (38) in the vicinity of the 
thresholds of the excitation of the tangling clustering in- 
stability. 

Let us consider the turbulent diffusion range of scales, 
l/\/Sc ^ i? <C 1 and introduce the following function 
^{t, R) = $(t, R) R^+^', where 



1 



11 = 



1 + 3(T^ 



1-2(T^ 



lOcTv 



1 



(49) 



Equation (38) is reduced to the Schrodinger type equa- 
tion: 



d^{t,R) 



-W^ + I{R), 



dt M{R) 
with the potential U and 1/M{R) in the form: 



(50) 



40crv 



1 

+(2cT^ 
1 



1 



5crv 



(1 + a^) 



1)1 



(1 + 

(2C7^ 



(Tv) (1 -I- 3(7^) 



(51) 



1 -|- 3(T.j, 

2 (l + 3a ) 2 7(ij) = Bi?i+M (52) 

M{R) 3 (1 + cr^) ' ^ ' ' 

and B = 20c7v/(l + CTv). Here we took into account 
that the main contribution to the source term I{R) 
for large Reynolds numbers is due to the first term 
B{R)N'^ in Eq. (23). Other contributions are negligible 
~ lORe-^/2 In^ Re [see Eqs. (15) and (37)-(35)]. 

We choose the initial conditions which correspond to 
a turbulence without particle clusters: '^{t = 0,R) = 



8 



0. Now we seek a solution of Eq. (50) in the follow- 
ing form: ^ YZ^x ! m{t)'^ m{R) and I{R) = 

B An^m(^), wlicrc „i{R) are the eigenfunctions 

determined by the following equation with I{R) = 0: 



1 



-*l-(^^ + 7m)*m = 0. 



M{R) ^"^ ^" ' ^^^^ 

Concequently, the function fm{t) is determined by the 
following equation: 



dt 



(54) 



Taking into account the orthogonality of the eigenfunc- 
tions, 



J^MiR)^m{R)^niR)dR 



(55) 



J,^M{R)<i>l{R)dR 

and solving Eq. (54) we obtain the expressions for the 
function fmit) and Am- 

20(7v A„ 



fm{t) = 



(1 + 0-v)7rr 



[exp(7„f) - 1] , 



_ J^M{R)R^+t^^miR)dR 
!^M{R)^1{R)dR ■ 



(56) 
(57) 



Taking into account 

[^m{t,R) R-i^+t^n 



that *™(i? 



1 we obtain the correla- 



tion function $(i? — ap,t) in the form: 
$(i? = ap,t) 



m=l ^"^ 



When the tangling clustering instability is excited, it 
causes formation of a cluster with the particle number 
density inside the cluster, which is much larger than the 
mean particle number density. 

The solution for $(t, i?) [sec Eq. (47)], which satisfies 
the above conditions has the following dimensional form: 



^t,R) = N'^ 



20fTv 



E 



Am 



cos 



(59) 



l + CTv 

X [exp(7mt) 

Here the correlation function $(i, i?) has the global max- 
imum at R = io/f-rp i-c, we assumed that in the molec- 
ular diffusion region the correlation function ^{t,R) is 
nearly constant, ^{t,R) « 1 [see Eq. (43)]. The first 
minimum of the correlation function <l>(i, R) for the mode 
m = 1 is located at i? = .Rmin = exp(l/A), and it is given 
by 



A^2 



= _:^Sc-V2 [exp(7it)-l] 
671 A 



where we took into account that cos [Kln(i?min/^z))] = 
—k/X. On the other hand, the maximum value of the 
correlation function ^{t,R) is 



N2 



BAi 
71 



[exp(7it) - 1] , 



(61) 
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FIG. 3. The particle number density inside the cluster 

nm3.^/N versus the particle radius Op for different values of 

n6. 



parameter F, and 
is given in /im. 



1, So = 10 Op. The particle radius Op 



where $n 



$(i? = ap,t) [see Eq. (58)]. This yields: 



eA I InSc 



(62) 



Since Up = N + n > 0, the function ^{t,R) = 
{n{t,x)n{t,y)) = {np{t,x)np{t,y)) - > -N^. There- 
fore, the minimal possible value of the function $(t, i?) 
is ^>min = — A^^. This condition together with Eq. (62) 
allow us to estimate the maximum number density of 
particles attained inside the cluster during the tangling 
clustering instability: 



N 



= |l + ^Sc 

TT 



A/2 



InSc 



1/2 



(63) 



The maximum value of the particle number density 

inside the cluster, riynax/N, versus the particle radius is 
shown in Fig. 3 for different values of parameter F. The 
discontinuity of the first derivative of rimax/N which is 
seen in Fig. 3 is related to the transition from one mech- 
anism of particle tangling clustering due to the source 
term to another mechanism caused by the tangling clus- 
tering instability. The exponential growth at the linear 
stage of the instability is saturated by the nonlinear ef- 
fects. The values of nmax/-^ in Fig. 3 are calculated using 
Eq. (63) that takes into account possible saturation of the 
tangling clustering instability caused by the exhaustion 
of the particles in the region surrounding the cluster. In- 
spection of Fig. 3 shows that the particle number density 
inside the cluster can increase in 10^ times compared to 
the mean particle number density. 

There are also other mechanisms leading to the nonlin- 
ear saturation of the tangling clustering instability dis- 
cussed in detail in Ref. 15. However, as follows from our 
analysis the main significant mechanism of saturation of 
the growth of the tangling clustering instability for small 
particles is exhaustion of the particles in the surrounding 
area. Indeed, the tangling clustering instability causes 
strong redistribution of particles so that inside the clus- 
ters the particle number density is strongly increased at 
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FIG. 4. The exponent A versus the particle radius ap for 
different values of parameter V and = 1. The particle 
radius Op is given in ^xa. 



FIG. 5. The particle number density inside the cluster 
fimax/N as a function of time for different values of the parti- 
cle radius Op and F = 10*, tr^ = 1, Sc = ICPup. The particle 
radius ap is given in /im. The dashed line is for > 3/ttm. 



some instant, while in surrounding regions it is essen- 
tially clccroasod. With the decrease of the number den- 
sity of particles the hydrodynamic description becomes 
inapplicable. It should be noted that we consider situ- 
ation when there is only the redistribution of the parti- 
cles without their creation or annihilation. This implies 
that particles from the cluster vicinity are concentrated 
in the central part of it, which can be expressed using the 
Corrsin integral of the correlation function of the parti- 
cle number density fluctuations: /p ^{t,R)R'^ dR = 0. 
This condition implies that the tail of the correlation 
function ^{t,R) must be negative [i.e., there is the anti- 
correlation tail of the function The transition 
from central positive part ^{t,R) > to the negative tail 
of <^{t,R) occurs at the distance that is of the order of 
several Kolmogorov scales. Note, that in contrast to the 
inertial clustering, the tangling clustering instability ac- 
cumulates particles into the cluster from the scales which 
are much larger than the Kolmogorov length scale. This 
is the reason why the concentration of particles inside the 
cluster increases due to the tangling clustering instability 
by several orders of magnitude. 

The value rimax depends strongly on the Schmidt num- 
ber Sc and on the exponent A [see Eq. (63)]. The expo- 
nent A depends on the degree of compressibility of the 
particle velocity field, cx St2r(Ma,Re,4/^T) see 
Eq. (35)]. The calculated values of exponent A versus 
the particle radius for different values of parameter F are 
shown in Fig. 4, that explains strong dependence of the 
maximum particle number density inside the cluster rimax 
on the particle radius and the parameter F. 

In general, other nonlinear mechanisms may limit the 
growth of the tangling clustering instability in the non- 
linear stage of its evolution, and therefore they limit 
the maximum attainable value of the particle number 
density, n'^'^^, attained inside the cluster. For exam- 
ple, a momentum coupling of particles and turbulent 
fluid becomes essential when the mass loading parameter 
<7 = mpTimax/p is of the order of unity. [1] Introducing 
a mean particle mass density Pp = nipN (e.g., in cloud 



physics it is a liquid water content measured in g/cm'^), 
we obtain the following constraint: n^^^/N < <,p/'Pp. 
Using p = 1.3 X 10^'^ g/cm'^ and Pp = 10~^ g/cm"^ we 
arrive at n^'^''/N < 1300 for c = 1. This limiting effect 
should be taken into account together with the satura- 
tion mechanism caused by the exhaustion of the particles 
in the surrounding area of the cluster. 

Figure 5 shows the temporal evolution of the number 
density of particles inside the cluster during the excita- 
tion of the tangling clustering instability for particles of 
different radius. Figure 5 reveals several interesting fea- 
tures pertinent to the instability, which deserve to be 
mentioned. The tangling clustering instability is less ef- 
fective for very small particles, ttp < 0.5/im. The reason 
is that this instability is saturated by the exhaustion of 
the particles in the surrounding area of the ehister at a 
very low value rimax/N ~ 13.2. For the paitieles hav- 
ing sub-micron and micron sizes the concentration inside 
the cluster can increase due to the tangling clustering 
instability up to several orders of magnitude. On the 
other hand, as follows from Fig. 2 the growth rate of the 
tangling clustering instability is the same for particles 
with flp = 0.565/im and for all particles with ap > 3/im. 
However, contrary to the case Up > 3/zm, the saturated 
value of the instability for ap = 0.565/im is very low, 
rim^/N « 27.2. 



VI. EFFECT OF DROPLET EVAPORATION ON 
TANGLING CLUSTERING INSTABILITY 

Let us study the effect of droplet evaporation on tan- 
gling clustering instability in stably stratified turbu- 
lence. The equation for the instantaneous number den- 
sity np{t,x) of droplets of the radius Up is given by: 

^ + V-(np v) = D„ Arip - ^ + /o, 

(64) 
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where the second term in the right hand part of Eq. (64) 
takes into account the droplet evaporation with the char- 
acteristic time Tev, and the last term, Iq, describes source 
of droplets due to condensation, which for simplicity is 
assumed to be constant. The equation for fluctuations of 
the droplet number density, n{t,x) = np(t, x) — N{t,x), 
reads: 

— + V-(nv- (nv)) =DmAn- {vV)N - NV-v 

(65) 



n 

Tev 



The last term in Eq. (65) describes the droplet evapora- 
tion. Using Eq. (65) we derive equation for the evolution 
of the two-point second-order correlation function of the 
droplet number density, $(t,R), see Eq. (17), in which 



/•OO 

D.^ (R) « 2 / {vi [0, 4(t, x|0)] vj [t, ^{t, X + R|r)] ) 
Jo 

xG{T)dT, (66) 

B(R)«2 / (6[0,^(t,x|0)]6[T,^(t,x + R|r)]) 
Jo 

xG{T)dT, (67) 

roc 

U,{K)^-2 / («,[0,C(t,x|0)]6[r,$(t,x + R|r)]) 

xG{T)dT, (68) 

G(t) = exp(— t/tcv), and other quantities are not 
changed by the droplet evaporation. 

Equation for $(<;, R) can be rewritten in a dimension- 
less form: 



9$ 

'dt 



M{R) 



B{R)- 



2r„ 



$ + J(i?), 



(69) 



where distance R is measured in units of Kolmogorov 
scale and time t is measured in units of /t^^ = 
Here the effective time t'^^ is 



(70) 



a modified turbulent diffusion time is = £'^/D^ , the 
effective turbulent diffusion coefficient in the Kol- 
mogorov scale is given by [50] 



and a modified function 1/M{R) is 



M{R) 



^ + ^-[l-F-{RF^y], 



(71) 



(72) 



where Sc*^^ = Scr'^^^/r^. Equation (69) shows that the 
evaporation decreases the source term [B{R) — 2T^/Tey\ 
of the generation of fiuctuations of the droplet number 
density. 
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FIG. 6. The growth rate 71 of the tangling clustering insta- 
bility (in units of 1/tj,) of the first mode (m = 1) versus the 
particle radius Op for different values of the relative humidity 
<p, and r = 10*, = 1, So = l(fap. The particle radius ap 
is given in fjcm. The dashed line corresponds to the inertial 
clustering instability. 



The analysis similar to that performed in Sect. 4 yields 
the growth rate for the mode m of the tangling clustering 
instability: 



1 



3(l-F3(T.j,) 
27r2m2(l 



200crv(cr^ - CTv) (3-Cr^)^ 



(l + av)2 



2(1 + C7^) 



2t2 



(73) 



where m = 1,2,3,.... This growth rate 7,,, of the sec- 
ond moment of particles number density was obtained 
by matching the correlation function ^{R) and its first 

derivative ^'{R) at the points R = l/VSc"^ and R= 1, 
which also yields: k/2(Ci -I- C2) « Trm /InSc^"'. There- 
fore, the evaporation of droplets causes decrease of the 
instability growth rate. Figure 6 shows the growth rate 
71 of the tangling clustering instability for different val- 
ues of the relative humidity (f> versus the particle radius. 
Here we used an expression for the evaporation time of 
droplets r'"' = 2.1 x 10~^ap/(l - </>), where the droplet 
radius is in microns and time is in seconds (see, e.g., 
Rcf. 51). As seen in Fig. 6, the evaporation of droplets 
drastically change the tangling clustering instability for 
small droplets, i.e., it increases the instability threshold 
in the droplet radius depending on the relative humidity 
(J). In addition, there is sharp maximum of the growth 
rate for 1 fim droplets if relative humidity is close to the 
supersaturated values: 99.8 % and 100 %. For low rela- 
tive humidity the growth rate of the tangling clustering 
instability is less in comparison with the supersaturated 
case. 

Figure 7 shows the temporal evolution of the number 
density of particles inside the cluster during the excita- 
tion of the tangling clustering instability for the relative 
humidity ^ = 90% and different droplet radius. This fig- 
ure also demonstrates that the evaporation of droplets 
drastically change the tangling clustering instability. 
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FIG. 7. The particle number density inside the cluster 
nmax/N as a function of time for different values of the par- 
ticle radius Op, and the relative humidity <f> = 90%, F = 10*, 
= 1, Sc = 10® Op. The particle radius Op is given in fj,m. 



VII. DISCUSSION AND CONCLUSIONS 

The present study has been inspired by the previous 
work [44] , where it was shown that the tangUng clustering 
of inertial particles in the temperature stratified turbu- 
lence holds the potential to promote a strong clustering 
with the considerably enhanced particle concentration in- 
side the cluster. Based on the thorough theoretical analy- 
sis, it is demonstrated in this study that the temperature 
fluctuations strongly contribute to the tangling cluster- 
ing instability. Temperature fluctuations caused by tan- 
gling of the mean temperature gradient by the velocity 
fluctuations, produce pressure fluctuations and enhance 
considerably particle clustering. The growth rate of the 



tangling clustering instability is in \/Re (Iq / Lt)"^ / {3Ma,)^ 
times larger than the growth rate of the inertial cluster- 
ing instability. 

The growth rate of the tangling clustering instability 
and the particle number density inside the cluster after 
saturation of the instability on the nonlinear stage de- 
pends on the parameter T w (Re^^^/dUa^^ 4 |VT| /T. 

We also found that depending on the parameters of tur- 
bulence and the mean temperature gradient there is the 
preferential clustering of particles of a particular size 
(the growth rate of the tangling clustering instability has 
a sharp maximum at this size). We demonstrated the 
strong effect of the droplet evaporation on this instabil- 
ity. The tangling clustering instability in the tempera- 
ture stratified turbulence may enhance significantly the 
collision rate of small particles, which is of interest for at- 
mospheric physics and many other practical applications. 
In particular this effect can substantially accelerate the 
coalescence of small droplets in atmospheric turbulence 
with temperature gradients. 
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